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The electromagnetic response of a pinned Abrikosov fluxoid is examined in the 
framework of the Bogoliubov-de Gennes formalism. The matrix elements and the 
selection rules for both the single photon (emission - absorption) and two photon 
(Raman scattering) processes are obtained. The results reveal striking asymmetries: 
light absorption by quasiparticle pair creation or single quasiparticle scattering can 
occur only if the handedness of the incident radiation is opposite to that of the 
vortex core states. We show how these effects will lead to nonreciprocal circular 
birefringence, and also predict structure in the frequency dependence of conductivity 
and in the differential cross section of the Raman scattering. 
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Recent far infrared transmission experiments on superconducting YBCO thin films [jT],|^ 
have revived once again the interest in the quasiparticle excitations of the vortex core. A 
few decades ago, the microscopic structure of a vortex hne in a type-II superconductor 
was the subject of many theoretical investigations |0-§1. Despite the rich structure in the 
density of states, predicted by the theory, early experiments on the microwave surface resis- 
tance ||2^ and microwave transmission of the mixed state confirmed only a very simple 



theoretical prediction: that the vortex core states occupy a fraction {H/Hc2Y of the total 
volume and that the vortex line is roughly equivalent to a normal cylinder of radius ~ ^, 
the superconducting coherence length . This was largely due to the lack of spatial and fre- 
quency resolution: the experimental signal was an average over a large number of vortices. 
More recently, however, rather spectacular results were obtained, by scanning tunneling 
microscopy (STM) ^^-^ which the spatial resolution of the signal was well below the 



coherence length. These experiments prompted theoretical work by various groups [[1^ -[T5 
, and it turned out, that one can understand very well the experimental results within the 
Bogoliubov - deGennes (BdG) formalism . 

Although the STM experiments reveal many aspects of the electronic structure of the 
vortex, they do not allow for measurements that probe the discreteness of the quasiparticle 
spectrum within the vortex, or the handedness imposed by the external magnetic field on the 
quasiparticle excitations of the vortex core. As we will show later on, this handedness will 



give rise to magneto-optical effects. It was pointed out to us by Girvin [0 and Kapitulnik 



r?! , that a high resolution optical gyroscope [0 might be able to detect the optical activity 
of the vortex core states. In order to probe the quasiparticle spectrum, Karrai et. al. 
have recently performed far infrared transition measurements on Y Ba2Cu^0j thin films in 
perpendicular magnetic fields. They observe a broadened edge-like feature at TTcm^^ and 
attribute this to the quasiparticle pair creation process inside the vortex core. Treating the 
mixed phase as a heterostructure, they give a phenomenological description of their data, 
and extract the resonant frequency mentioned above, which, however, leads to a somewhat 
large energy gap : A = 63meV. 



In this paper we examine the response of a pinned Abrikosov vortex to an external 
electromagnetic field within the BdG formalism. We are interested in the case of an isolated, 
static vortex line carrying a single fiux quantum. The external magnetic field is taken in the 
+z direction, with Hd < H <^ Hc2- Due to the cylindrical symmetry, one can choose the 
quasiparticle amplitudes in the following form { r = {t, 6, z )) B: 
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Here cTj 's are the usual Pauli matrices. The magnetic quantum number fi is half an odd 
integer. If we measure the energy relative to the fermi energy, Ep, and the length in units 
of ^, in a gauge where the pair potential A{x) is real the BdG equations read 
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In the above equation g is a two component spinor, C = h'^ / (2m^^) and k"^ = kp — k"^. (For 
simplicity we have assumed an isotropic effective mass.) In general, a self-consistent solution 
to the above equation can be obtained only numerically ||14| , p!5| . The general features of the 
solution are the following : the bound states, with exponentially decaying quasiparticle 
amplitudes, have an energy spectrum with E{fi) oc /i at small fi, the spacing between the 
levels being ~ A'^/Ep where Aqo is gap in the bulk. The negative energy states are fully 
occupied at T = 0, whereas the positive energy states are empty 0. The scattering states 
have a continuous spectrum, with energies \E\ > A^q. 

The perturbing hamiltonian describing single-photon (emission - absorption) processes, 
can be given as 
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where the sum is over the spin indices. The power absorbed by the system irradiated with an 
external electromagnetic field can be given in two different ways: V = (2^7^/0^)0-1(0;) = huW 
where <Ji{u!) is the real part of the conductivity <j{uj) = (ri{u)+i(T2{u!) and W is the transition 
rate of the system under the influence of the perturbation Tii. 



O-TT 

^ = T2 E \Hfr\^5{Ef -E,- hu) (4) 

where Hji = {f\T-Ci\i). Depending on the final state |/), the matrix elements will involve 
different coherence factors. For example, when the photon emission - absorption process 
scatters a single particle from one state to another the transition rate is given as 



We^.abs = ^ ( — 1 E Wn,n'? >< fn' (1 " f n) ^{En ' ± hu) , (5) 



where fn = f{En) is the Fermi function. We consider here the absorption process only, since 
we are ultimately interested in the low temperature limit, where there is no spontaneous 
photon emission. The most interesting case for the present problem is that of the circular 
polarized light : A{r) = Aqe± exp{i{qf ~ ut)) where e± stands for the usual polarization 
vectors corresponding to the left / right circular polarized light, e± = ±iey (note that we 
consider the case of a transverse electromagnetic wave propagating parallel to, and in the 
same direction as, the external magnetic field). In this particular case the matrix elements 
are 



— I dr 



dr ' 2' 

' } , (6) 



implying the following selection rules: 



k,-k'^-q = Q (7) 
/i' - ± 1 = (8) 
E„ -E^,-nuj = Q. (9) 

Further simplification occurs when we take the low temperature limit: then all the states 
with negative angular momentum /i' < are occupied and therefore the selection rules 
can be satisfied only for /i = /i' + 1 with the final state yU = +1/2 and initial state ^ = —1/2. 
But this is possible only when the light is left circularly polarized . Explicitly inserting these 
values for the angular momenta, we find 



Kn' = / drr{g:;^^g:, - g-^^g',} (10) 

and 

M" , = . (11) 

Thus, the chirahty of the vortex core states becomes manifest in the above selection rules 
governing the absorption of circularly polarized electromagnetic radiation. 

The pair creation process, where the incoming photon creates a pair of quasiparticle 
excitations, can be discussed in a similar fashion. At zero temperature only the /i > are 
available, so that a photon with an energy just above the pair creation threshold will create 
a pair of quasiparticles with opposite spins but the same angular momentum /i = 1/2. The 
transition probability corresponding to this process can be obtained by explicit calculation 
in the same way as that of the single particle process. It is, however, possible to obtain 
the selection rules and the matrix elements corresponding to this process directly from the 
single particle quantities, via a particle - hole like symmetry connecting the positive and 
negative energy states ||20|. According to this symmetry |Q, if ((7^,(7^) describes the state 



with energy and momentum En-, fcz, /i, then {g~, —g^) corresponds to —En, —k^, —^i- In this 
way, the single particle process, where a negative energy state is destroyed and a positive 
energy state is created, can be easily related to the pair creation where two positive energy 
states are created. In particular, the selection rule for the angular momentum for this process 
is that fi' + fi = ±1. Since the pair with the lowest possible energy has /i' = /i = 1/2, the 
selection rule can again be satisfied only by a light carrying positive helicity, i.e., having left 
circular polarization. 

Similar results were obtained independently, and parallel with our work, by Zhu, Zhang 
and Drew [^. These authors consider the superconducting ground state as a vacuum. In 



this picture, which is equivalent to ours, the electromagnetic radiation can be absorbed only 
by pair creation at T = 0. They also point out that if the carriers in a type II superconductor 
are hole type, such as in most of the high Tc compounds, the situation is reversed, with only 



right circularly polarized light is absorbed at zero temperature. This is a consequence of the 
CT invariance (simultaneous charge conjugation and time reversal). 

These asymmetries will probably lead to experimentally observable consequences for the 
following reason. A difference in conductivity will cause a difference in the refractive index 
of the system with respect to the light waves of different circular polarization, since the 
complex refractive index is given by N'^{uj) = Coo + i4:TTcr{uj)/uj = {n + m)^ where eoo is the 
dielectric constant at large frequencies, n is the real refractive index and k is the absorption 
coefficient. Different refractive indices for the two circular polarizations will result in 
a nonzero Faraday angle, with which the polarization plane of a linearly polarized light is 
rotated by a sample with thickness z : (pp = {(jjz/2c){n^ — nJ). 

An attempt was made by Karrai et. al. p[ to check the angular momentum selection 
rule and look for chirality in the vortex response by performing circularly polarized light 
transmission measurements on superconducting Y Ba2Cu^0'j thin films. They found no 
evidence for optical activity in the vortex response. Instead, the signal is dominated by 
magneto-optical effects attributed to the condensate and it is interpreted as a cyclotron 
resonance of the superconducting ground state. 

A possible explanation for the lack of optical activity has been proposed by Hsu |T^. In 



this alternative picture, the /i = — 1/2 ^ +1/2 dipole transition is hidden by the resonance 
of the circular motion performed by an unpinned vortex, which occurs at the same frequency 
as the dipole transition ~ A^^/Ef. This resonance is clearly not sensitive to the polarization 
of the external electromagnetic radiation, which drives the motion of the vortex. However, 
we find it surprising that the resonance of this circular motion would occur at such a high 
frequency. 

Let us now discuss the inelastic light scattering on the vortex core states. The electronic 
Raman scattering in metals with energy and momentum transfer uj = uJi—uJs and q = ki — kg 
can be understood as a scattering on the effective density pg = J2k,a'lkak+q,a^k,a p2|j23 



where the scattering strength 7^ is strongly polarization dependent and satisifies 7^ = 7_fc 
by time reversal symmetry. Note that the continuity equation is not necessarily valid for pq 
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2J]. In direct space this is equivalent to the following effective density: 

Pe//(r) = / d'RJ2^|Ji{r)Mr + Rh{R) • (12) 

For an isotropic system 7^ = const and consequently 7(-R) = S{R). The photon cross section 
per unit area and time is given by p5| : 



5^ - h (S)' POT^ J:\im.nm - E. - n.) . (13) 

where 6 is the angle of incidence and the effective interaction hamiltonian is defined as 

^e// = "^AkAk, [ d^rdV e"^^{r - r)^t (^)^^(/) . (14) 

Here tq = e^/mc^ is the Thomson radius. We are now able to examine the matrix elements 
corresponding to different processes. For example, a photon can be absorbed and re-emitted 
by a particle, or a pair can be created which will finally recombine into the condensate and 
provide the outgoing photon by emission. The fermion subspace of the initial and final state 
is practically the same as for the single photon processes we've discussed. In the case of 
an isotropic system, with 7(r) = 6{r) the matrix elements can be easily calculated. Here, 
for the sake of simplicity, we choose a special experimental geometry in which the direction 
of the momentum transfer q is along the z axis. The generalization for other geometries is 
quite straightforward. The single particle processes are described by: 

Mn,n' = J drr{g+g+, - g-g-,} , (15) 
along with the selection rules: 

h = k', + q (16) 

fi' = fi (17) 

= E^i + hu . (18) 

The corresponding relations for the pair creation process can be obtained with the particle 
- hole symmetry mentioned before. 



The anisotropic case can also be discussed if we assume that the function 7(r) obeys 
cyhndrical symmetry. This is approximately the case in a layered superconductor, when 
the axis of the vortices and the incident and the scattered light waves are perpendicular 
to the layers. Then one can choose the following trial form : 7(r*) = 7m,r('")e*'"^e*'"^. The 
evaluation of the matrix elements in this case is somewhat lengthy but still straightforward. 
As a result we obtain that the selection rules are still the same as in the isotropic case : this 
is a consequence of the cylindrically symmetric scattering strength 7(r). 

The more interesting feature is the dependence of the scattering cross section on the 
transfer energy. This will be investigated by evaluating the matrix elements explicitly, 
using the numerical solutions [Il4| , pr5| of the Bogoliubov - de Gennes equations. The detailed 
numerical results will be reported elsewhere. Nevertheless, it is possible to discuss intuitively 
the most important features of the frequency dependence. In a bulk superconducting phase 
the dynamic structure factor is zero until the energy transfer reaches the huo = 2 A threshold : 
the superconducting ground state cannot absorb energy unless the photon has enough energy 
to break a Cooper pair. Above the threshold, the dynamic structure factor is divergent: the 
states excluded from the gap pile up above 2A giving rise to a square root singularity in 
the density of states and consequently in the dynamic structure factor. This divergency 
is in practice removed by gap anisotropy or final state interaction In a type 



II superconductor in the mixed phase, the situation is quite different: the quasiparticle 
excitations of the vortex core are capable of absorbing the energy of the incoming photon. 
As mentioned before, these states occupy a fraction {H/ i/c2)^ of the total volume and provide 
a quasi-normal region where the dissipation is possible unless the energy of the incoming 
photon is less then twice the minigap : huj < 2A^/Ef. Although this minigap is often 
very small in conventional superconductors, it is larger in a few materials. For example, in 
NbaSn, 2A'^ / Ep ~ 0.97^ which is is not necessarily unattainable experimentally p7| , |28| . For 
high Tc superconductors, this value is even higher. 

It is also interesting to ask whether the interaction between the quasiparticles would 
have significant effect on the Raman spectrum. If the final state interaction is included in 
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the description of Raman scattering on superfluid He II or superconductors [^] , the 



quasiparticles can form a bound state below the threshold. Similar effects may occur in the 
mixed state of superconductors. 

In conclusion, we have investigated the electromagnetic response of the superconducting 
vortex states. Due to the handedness of the vortex states, strong asymmetries are discovered 
for the low temperature absorption of circularly polarized light : the fi = —1/2 —* fi = 
1/2 dipole transition and quasiparticle pair creation are possible only with left circular 
polarization, a striking consequence of the handedness of the vortex core states. We related 
the matrix elements to experimentally accessible quantities such as the dissipative part of 
the conductivity and the Faraday angle, and argued that the asymmetry in the absorption 
of circularly polarized light in a superconducting vortex and a finite rotation of the plane 
of polarization of a linearly polarized light at low temperatures might be experimentally 
observable. 

We proposed that inelastic light (Raman) scattering could be used to investigate the 
vortex core states and developed the theory of Raman scattering on a superconducting 
vortex for isotropic and two dimensional, layered superconductors. We argued that the 
presence of the vortex core states will lead to finite Raman intensity below the usual 2A 
threshold and that the minigap in the discrete spectrum of the vortex core states may be 
observable using low frequency Raman spectroscopy. This is even more likely in high Tc 
superconductors, where such a measurement, if successful, can reveal the magnitude of the 
energy gap. The results of conventional measurements for the energy gap in the cuprates, 
such as tunneling, infrared and high frequency Raman spectroscopy are not yet consistent, 
mainly due to the fact that at high frequencies the spectrum is convolved with the phonon 
and other excitation spectra [ pO|] . 
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